Unlike Lennard- Jones Parameters for Vapor-Liquid Equilibria 

Thorsten Schnabel, Jadran VrabecQ, Hans Hasse 
Institut fiir Technische Thermodynamik und Thermische Verfahrenstechnik, 
Universitat Stuttgart, D-70550 Stuttgart, Germany 

Abstract 

The influence of the unUke Lennard- Jones (LJ) parameters on vapor-Hquid equiUbria of 
mixtures is investigated and the performance of eleven combining rules is assessed. In 
the first part of the work, the influence of the unlike LJ size and energy parameter on 
vapor pressure, bubble density and dew point composition is systematically studied for the 
mixtures CO+C2H6 and N2+C3H6, respectively. It is found that mixture vapor pressure 
depends strongly both on the size and the energy parameter whereas the bubble density 
depends mostly on the size parameter and the dew point composition is rather insensitive 
to both parameters. In preceding work, unlike LJ parameters were adjusted to experimental 
binary vapor-liquid equilibria for 44 real mixtures. On the basis of these results, in the 
second part of the work eleven combining rules are assessed regarding their predictive 
power A comparison with the adjusted unlike LJ parameters determined from the fit shows 
that none of the eleven combining rules yields appropriate parameters in general. To obtain 
an accurate mixture model, the unlike dispersive interaction should therefore be adjusted to 
experimental binary data. The results from the present work indicate that it is sufficient to 
use the Lorenz rule for the unlike LJ size parameter and to fit the unlike LJ energy parameter 
to the vapor pressure. 
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1 Introduction 



In molecular simulations of a binary mixture A+B with pairwise additive potentials, three dif- 
ferent interactions occur: two like interactions between molecules of the same type A-A and 
B-B, which are fully defined by the pure component models, and the unlike interaction between 
molecules of different type A-B. In mixtures consisting of polar molecules, the electrostatic part 
of the unlike interaction is fully determined by the laws of electrostatics. However, there is no 
rigorous physical framework that yields reliable unlike repulsion and dispersion parameters like 
the Lennard- Jones (LJ) parameters studied in the present work. For finding these parameters, 
combining rules were developed in the past based on physical and mathematical intuition or on 
empirical approaches. Eleven of these combining rules were investigated in the present work. 
These combining rules rely soley on pure component data, namely the LJ parameters and, in 
some cases, additionally the polarizablility a or the ionization potential /. Other combining 
rules, that are not discussed in this work, also employ dispersion force coefficients |lll|2l[3]|, 
diamagnetic susceptibility ijHO or effective transition energies (61171. 

Another approach for obtaining unlike LJ parameters is to adjust them directly to experi- 
mental binary data, for which a single data point may in principle be sufficient. Kohler et al. dH 
fitted both unlike LJ parameters to experimental second virial coefficients of binary mixtures. 
Following this approach, MoUer et al. developed a method for adjusting the unlike LJ size 
and energy parameters to experimental excess volumes and enthalpies. Their investigations 
showed that the unlike LJ size parameter determined from the fits practically does not deviate 
from the arithmetic mean of the pure component LJ size parameters. This finding is supported 
by the work of Vrabec and Fischer for Ar + CH4 |jTOl as well as for the three binary mixtures 
that can be formed out of CH4, C2H6 and CO2 [HJ- Kronome et al. [fT2ll applied a similar 
approach, specifying the unlike LJ size parameter by the arithmetic mean. They adjusted only 
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the LJ energy parameter to the experimental excess Gibbs enthalpy for the mixture N2+C2H6 
and obtained favorable results for vapor-liquid equilibria. 

In previous work of our group |fT3l[T4l[T5l . unlike LJ energy parameters of 44 binary mix- 
tures were fitted to one experimental binary vapor pressure each. These mixtures contain noble 
gases, homonuclear and heteronuclear diatomics, small hydrocarbons, carbon dioxide, carbon 
disulfide and halogenated hydrocarbons, i.e. refrigerants. Overall 22 components are studied, 
where the molecules are composed out of 1 to 9 atoms. These components were modeled by 
the one-center Lennard- Jones potential (ICLJ) |fT6l , the symmetric two-center Lennard- Jones 
potential either with a pointdipole (2CLJD) IfTTll or with a linear elongated pointquadrupole 
(2CLJQ) [16]. The parameters of the pure component molecular models were adjusted to ex- 
perimental pure component bubble density and vapor pressure data |fT6l [TTl . The pure com- 
ponent models are accurate, the mean errors of the vapor pressure, bubble density and heat of 
vaporization are typically 4 %, 1 % and 3 %, respecitvely. Vapor- liquid equilibria of the 44 mix- 
tures are described with typical deviations in vapor pressure, bubble density and in dew point 
composition of below 5 %, 1 % and 0.02 mol/mol, respectively lfT3l [T4l [TSl . These mixture 
models also predict favorably other fluid properties like Joule-Thomson inversion curves [18] 
or self-diffusion and binary Maxwell-Stefan diffusion coefficients as well as bulk viscosities 
and thermal conductivities |[T9ll20ll . 

Ternary mixtures are predicted without further parameterization since the force fields are 
based on pairwise additive interactions. Results for ternary vapor-liquid equilibra of six mix- 
tures [[T3l [Ml [151 |2T]| confirm the predictive power of such models. 
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2 Combining rules 



The LJ potential u}^^ is the most widely used functional form for describing repulsion and dis- 
persive attraction. This pairwise additive potential acts between two molecules / and j and is 
given by 
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where a is the site index of molecule i and h the site index of molecule j. ml and ml are the 
number of LJ interaction sites of molecule i and j, respectively. The site-site distance is rijab^ 
Oab and Cab are the LJ size and energy parameters between sites a and b. For obtaining the LJ 
parameters of the unlike interaction, denoted by Oab and E^^ subsequently, usually combining 
rules are used. 

Lorentz-Berthelot (LB) 

The Lorentz-Berthelot combining rule is most widely used to determine the unlike LJ parame- 
ters. Lorentz f2l\ proposed to use the arithmetic mean for the unlike size parameter motivated 
by collisions of hard spheres 



LB _ '^aa + '^bb 
^ab — ^ • 



Berthelot 11231 proposed with little physical argument the geometric mean for the unlike 
energy parameter 



Eflf = V^aa^bb- (3) 

The LB combining rule is by far the oldest and most common approach, but due to the fact 
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that it can lead to inaccurate mixture properties [|24ll25l . numerous other combining rules have 
been developed. 

Kohler (K) 

Kohler [|26]| used the approach of London Il27ll28]| for the dispersion energy to derive a combin- 
ing rule for the unlike LJ energy parameter which uses the polarizability a 



^ab — ^ 



<3aa<3bb 



^aa^bb ■ (4) 



^Ib^aa^aa + al^oli^ehb 



Oaa+^bb. 

The application of polarizability in combining rules is not straightforward. Generally, the 
polarizability quantifies the distortion of the overall charge distribution of an atom, group or 
molecule by an electric field. It is in fact a tensor, which includes the anisotropic contributions 
of the directionally distorted charge distribution. This tensor is often described by a scalar a 
which is appropriate to characterize the distortion sensitivity for the whole molecule. Laidig 
and Bader ||29l suggest that atomic or group polarizabilities influenced by the chemical enviro- 
ment additively contribute to the overall molecular polarizability. Hence, the atomic or group 
polarizabilities required to apply the combining rules to the symmetric two-center LJ potential 
ascribe 50 % of the molecular polarizability to each of the two LJ sites. For molecules modeled 
with the one-center LJ potential, e.g. methane, the total molecular polarizability was used here. 
Polarizabilities were taken from [|30l except for the refrigerants R125 |l3ll|32l, R134a [[SllEIl 
andR152a|l32l[3l. 

Furthermore, Kohler proposed to use the arithmetic mean for the unlike LJ size parameter 
(Lorentz rule) as defined by Equation 

Hudson-McCoubrey (HMC) 

Hudson and McCoubrey [l34l reformulated Reed's proposal [|35l for the calculation of unlike 
LJ energy parameters which considers the ionization potential / of the pure components and 



volume effects of the mixture. They claim that their expression is "suitable to obtain viral 
coefficients of binary gas mixtures and binary gas-liquid critical temperatures" ||34l 
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Oaa + (3bb 

The ionization potential is defined as the energy which is necessary to remove the outermost 
electron from an atom or molecule. Since the outermost electron cannot be explicitly assigned 
to a specific atom or group of a molecule, the full molecular ionization potential is ascribed in 
the present work to every LJ site of the molecular model. Ionization potentials were taken from 
[|30l except for the refrigerants R125 Bill [21, R134a and R152a [|32ll33l. 

Srivastava and Madan [|36ll derived the same expression as Equation © with an approach 
similar to that of Kohler [|26l . Note that Equation dH) can be transformed to Equation ([5]) if the 
following relation between polarizability, ionization potential, LJ size and energy parameters is 
assumed for both components a and b 



^^aa^aa ~ ^ ^aa^aa ■ (6) 

Some authors, e.g. Good and Hope Q or Reid and Leland [l37ll . suggest that the ratio between 
the geometric and arithmetic mean of the ionization potentials in Equation Q can be set to 
unity. However, in the present work Equation ^ was applied. 

For the unlike LJ size parameter in the HMC combining rule the Lorentz rule is applied as 
defined by Equation Q. 

Fender-Halsey (FH) 

Fender and Halsey [|38l proposed to use the harmonic mean for the unlike LJ energy parameter 
in their empirical combining rule 
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(7) 



which always gives smaller or equal values than the geometric mean. Mixtures of noble gases 
have mostly smaller unlike LJ energy parameters than those obtained with the geometric mean. 
Hence, Fender and Halsey [[38l achieved with their harmonic mean combining rule more accu- 
rate second virial coefficents for the mixture argon+krypton than with the LB combining rule. 

For the unlike LJ size parameter in the FH combining rule the Lorentz rule is applied as 
defined by Equation Q. 



Hiza et al. [|39l l40l 1411 used a semi-empirical approach to correct the LB combing rule. The 
expressions for the correction factors, dependent on the ionization potentials and empirical 
factors, were determined by fitting to equilibrium, transport and beam scattering data of noble 
gases and light hydrocarbons 
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and 
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with 
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where the component order has to be chosen so that laa > hb- 
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Since kah from Equation (flOl) is always larger than zero for nonidentical ionization poten- 
tials, the resulting unlike LJ size parameters are always larger than those from the Lorentz rule. 
The opposite holds for the unlike LJ energy parameters when compared to the Berthelot rule. 

Sikora (S) 

In Sikora's combining rule [|42ll repulsion is considered via deformation energies of the electron 
clouds and resulting unsymmetric collision diameters which occur at small distances between 
overlapping atoms or molecules 



Kt = 2 
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Applying similar considerations for the unlike LJ size parameter [|43l . the following equa- 
tion was proposed 
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Smith-Kong (SK) 

Smith [44J developed a combining rule considering the repulsive interaction of simple molecules. 
His combining rule includes the atomic distortion theory for repulsion. Kong [|45l used Smith's 
approach and a geometric mean relationship for the attractive interactions and applied it to the 
LJ potential leading to 
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Jeaaal2)T7+(£,,ol2)- 

Halgren (HHG) 

From investigations of noble gas mixtures using ab initio and experimental data, Halgren H6ll 
deduced that the experimental pattern of the van der Waals minimum energy distances is well 
described by the "cubic-mean" combining rule 



^HHG ^flfl + 
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^la+<^lb 



(15) 



The van der Waals minimum energy distance is associated with the LJ size parameter. To 
introduce the experimental unlike noble gas values into the LJ energy parameters, a combination 
of the geometric and harmonic mean was applied which is claimed to yield good descriptions 
of the experimental noble gas values 



„HHG _ 
^ab — 



^^aa^bb 



,2 ■ 



(16) 



Waldman-Hagler (WH) 

Waldman and Hagler HTl derived their combining rule by applying mathematical methods 
which consider symmetry, uniformal scaling and simplification of a two-parameter to a single- 
parameter problem to deduce the general functional form of their combining rule. From a 
"graphical analysis" of noble gas mixture data they obtained 
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^aa+< 



(17) 



and 
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Al-Matar (Ml and M2) 

Al-Matar and Rockstraw [|48l proposed combining rules which are based on similar mathemat- 
ical considerations as used by Waldman-Hagler (see above). The first of their combining rule, 
abbreviated in the following with Ml, was obtained by functional analysis of experimental no- 
ble gas mixture data. They consider the pure component LJ parameters equally weighted. Ml 
is given by 



^ L=0 \ ^aa^bb 



and 



^ab — —^y^aa^bb- (^U) 

2 /0.25(aL+aL)'^'"^" 



L 

L=0 V '-'"""bh 

The second combining rule, abbreviated in the following with M2, uses weighting matrices 
which can account for uneven contributions of the pure component LJ parameters to the unlike 
quantities. The weighting matrices of M2 were obtained by fitting them to noble gas mixture 
data. The unlike parameters are given by 

o^^ = {0.2S20al + 0A132alol + 0.24ASal) ' , (21) 

and 

^M2 _ 0.03995e^, + 0.9564698845v/eA + 0-00355efc^ ^3 ^3 
0.2820a6, + 0.4732o3„o3, + 0.2448o6, 



10 



In Equation (|2TI) . the lower valued LJ parameter has to be taken for site index a [|48l . since 
unequal weights are used. Table [U summarizes the eleven studied combining rules. 



3 Case Study 

Before studying the performance of the different combining rules, it is useful to get a picture of 
the influence of Oab and tab on the vapor- liquid equilibria of mixtures. Vapor pressure, bubble 
density and dew point compositon at specified temperature and bubble point composition were 
investigated regarding their sensitivity on (5ab and Zab- Instead of using the absolute values of 
Oab and Eab^ the unlike LJ parameters are discussed here in terms of the deviations T] and ^ from 
the Lorentz-Berthelot combining rule 



(23) 



and 



Eflfc = "^y/^aa^bb- (24) 

The influence of the unlike LJ parameters on the vapor-liquid equilibria of mixtures was 
investigated in a case study for which two different binary mixtures (N2-1-C3H6 and CO-1-C2H6) 
were chosen here. Three of the regarded components were described by 2CLJQ models, only 
for CO a 2CLJD model was used. For CO-1-C2H6, the case study was performed at 223 K and 
xco = 0. 198 mol/mol, where experimental vapor pressure and dew point carbon monoxide mole 
fraction are 5.614 MPa and 0.8065 mol/mol [49]. To eliminate the influence of experimental 
scatter, the pressure - composition data were smoothed [|50ll . The reference bubble density of 
16.499 mol/1 was determined using DDMIX [[5T1l which was designed to yield reliable mixture 
densities. 
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The experimental data for N2+C3H6 at 290 K and jcn2 = 0.1 146 mol/mol, where the case 
study was performed, are for the vapor pressure and dew point nitrogen mole fraction 6.783 MPa 
and 0.7497 mol/mol ||52]| . Again smoothed values were taken here IISOl . The bubble density of 
12.4094 mol/1 was estimated by the Rackett model lEll |54l together with the mixing rule of 
Chueh and Prausnitz ll55l . 

Vapor-liquid equilibria for CO+C2H6 were simulated for every combination of 
ri G [0.96, 0.98, 1, 1.02, 1.04] and ^ e [0.96, 0.98, 1, 1.02, 1.04], i.e. overall 25 combinations, 
using the Grand Equilibrium method [[56ll . Simulation results of vapor pressure, bubble density 
and dew point composition are therefore functions of T] and ^. Their deviations to the exper- 
imental values are illustrated as three-dimensional plots in Figures [H to [3] where the surfaces 
were derived from the 25 simulation points. All combinations of r\ and ^ which yield the exper- 
imental value are plotted as a solid slim line in Figures [T] to [3] and as a projection onto the T]-^ 
plane. For the further discussion experimental and simulation uncertainties are needed, where 
the following estimates were used: 3 % for the vapor pressure, 0.5 % for the bubble density 
and 4 % (about 0.03 mol/mol) for the dew point composition. These limits are also included in 
Figures [H to [3] as solid bold lines. 

The vapor pressure significantly dependens on both r\ and ^, cf. Figure [IJ It varies by about 
35 % for r\ (at constant ^) and 45 % for ^ (at contant r|) in the investigated parameter range. Due 
to that high sensitivity, the range where the simulated and experimental vapor pressure agree 
within 3 % is rather small. 

In contrast, the bubble density is practially independent of the parameter ^ as can bee seen 
from Figure [21 but it varies by about 6 % in the investigated range oft]. The experimental value 
can be obtained by choosing r| close to unity in combination with any value of ^ in the regarded 
range. 
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The dew point carbon monoxide mole fraction depends both on r\ and ^ as shown in Figure 
[3l But it varies only by about 6 % for r\ (at constant ^) and 2 % for ^ (at contant Tj) in the 
investigated parameter range. 

Summarizing Figures [T] to [3] in Figure IH it can be seen that there is one combination of r\ 
and ^ where both vapor pressure and bubble density of molecular mixture model and experiment 
coincide exactly (r\ = 0.9972 and ^ = 1.0145). With the uncertainties for the bubble density 
and the vapor pressure, a target area is defined in which the combination of r\ and ^ yields an 
accurate description of these both properties, cf. the shaded area in Figure IH The experimental 
dew point composition within its (large) uncertainty is also met in this target area. 

Figure |4] also allows an easy assessment of the combining rules introduced in Section [2l 
The LB combining rule is the closest to the target area, but outside of it. All combining rules 
shown, except SK, describe the bubble density and dew point composition within their assigned 
uncertainty. However, the experimental vapor pressure is predicted by none of them. The 
LB combining rule (tj = 1 and ^ = 1) yields deviations from the experimental vapor pressure, 
bubble density and dew point carbon monoxide mole fraction of +5.6 %, —0.3 % and +1.8 
%, respectively. The SK combining rule gives the most remote combination (r] = 1.0059 and 
^ = 0.9630) with deviations of +21.7 %, -0.8 % and +2.3 %, respectively. 

An analogous systematic study for N2+C3H6 includes 30 vapor-liquid equilibria simulation 
results for Tj G [0.96, 0.98, 1, 1.02, 1.04] and ^ G [0.94,0.96, 0.98, 1, 1.02, 1.04] at 290 K and 
xm = 0.1146 mol/mol. Since the experimental bubble density was estimated with the Racket 
model, a higher uncertainty of 1 % was assigned to it. Figure \5\ shows the comparison to 
experimental data where the target area is obtained analogously to the study of CO+C2H6. 
Only the M2 combining rule is within this area, whereas HMC is very close to it. For M2 
(T] — 1.0026, ^ = 0.9538) deviations from the experimental vapor pressure, bubble density and 
dew point nitrogen mole fraction of — 1.7 %, —0.8 % and +3.0 % were found, respectively. SK 
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(ri = 1.013, ^ = 0.9101) yields again the most remote results with deviations of +15.5 %, —1.6 
% and +4.8 %, respectively. 

From both case studies, the conclusion is drawn that the vapor pressure, which significantly 
depends on both unlike LJ parameters, is the most difficult property to be predicted by combin- 
ing rules. The bubble density is mainly determined by the unlike LJ size parameter and values 
close to the arithmetic mean (rj = 1) describe the experimental bubble densities well (see also 
[Sol [HI). This is achieved by most of the combining rules. The dew point composition depends 
on both r| and ^ like the vapor pressure, however, its sensitivity is considerably weaker com- 
pared to the vapor pressure. Due to these findings, it can be recommended to use the Lorentz 
rule for the size parameter since it predicts bubble densities well and to adjust the unlike energy 
parameter to the vapor pressure. 

4 Comprehensive Study 

In previous work ffT3l [T4l [T5l . molecular models for 44 binary mixtures were developed. For 
describing the mixture, the Lorentz rule was used for Oab, while Eab was determined using the 
adjustment parameter ^ as defined by Equation (|24l) . The state independent parameter ^ was 
optimized by a fit to one experimental vapor pressure data point of the mixture. That procedure 
chosen in [[T3l[T4l[T5l is impressively supported by the results of the case study presented above. 

The results from [fT3l[T4l[T5l can be directly used to assess combining rules. For combining 
rules which employ the Lorentz rule for Oat this is completely straight-forward as the adjusted 
^ parameter is known and only needs to be compared to the results for the different combining 
rules. But also combining rules that do not employ the Lorentz rule can be included in the 
assessment as can be seen from the results of the case study and explained below in more detail. 

The pure components of the 44 binary mixtures are of three different molecular model types: 
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spherical and non-polar (ICL J), anisotropic and dipolar (2CLJD) and anisotropic and quadrupo- 
lar (2CLJQ). Hence, six mixture types can be distinguished as given by Table [2l The adjusted 
unlike LJ energy parameters of those 44 mixture models deviate by up to ±10 % from the 
Berthelot rule, cf. Table [3] or see lfT3l[T4l[T5l. 

Seven of the eleven combining rules do not use the Lorentz rule for the unlike LJ size 
parameter. For an assessment, it has to be investigated by how much they deviate from the 
Lorentz rule. This was done here with the root mean squares (RMS) 



where A'^ is the number of mixture models to which the particular combining rules (x = H, S, 
SK, HHG, WH, Ml, M2) were applied. In Table H RMS^ are given distinguishing between the 
mixture types 1 to 6. It additionally incorporates data taking all 44 binary mixture models into 
account. 

Table |4] shows that the S combining rule yields the same unlike LJ size parameters as the 
Lorentz rule for all mixture types. Also M2 and SK are very close, deviating by approximately 
0.5 % only. The remainder, H, HHG, WH and Ml, yield all mean deviations in a narrow band 
at 1.4% and below. For mixture type 6, all seven combining rules yield results very close to 
the arithmetic mean due to fact that all components have very similar parameters. Figures |4] 
and [5] give an impression of the influence of a variation of T] of 1 % on the mixture vapor- liquid 
equilibria. For the pressure and the dew point composition that variation will typically result in 
changes that do not exceed the experimental uncertainty. For the bubble density the changes are 
more significant but still the changes do not exceed 1 % in the bubble density. Thus, it can be 
concluded that there is only a small influence of the combining rules via Oab on the vapor-liquid 
equilibrium. Therefore, in the following only results of the combining rules for ^ are discussed. 




(25) 
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Unlike LJ energy parameters predicted by the eleven combining rules are compared to the 
adjusted values in Figure [6l The adequacy of the investigated combining rules is indicated by 
the distance from the reference line (unity), which represents the adjusted LJ energy parameter. 
Generally, too low dispersion energies are predicted, which is particlularly visible for mixture 
types 2 to 5. Note that it is implicitly assumed in Figure [6] that the adjusted parameter for Eab 
determined in |fT3l [T4l [T5l using the Lorentz rule for Oab also holds for the combining rules 
which apply other approaches for determining Oab- 

To assess the predictive power of the eleven combining rules the deviations of the unlike LJ 
energy parameters of the particular combining rules from the adjusted values were summarized 
by the root mean square 



The deviations for all eleven combining rules, distinguishing between mixture types, are 
given in Table [51 For mixture type 1, for which most of the eleven combining rules were tested, 
they are better than the LB combining rule. Many of them were fitted to noble gas mixture 
data and they indeed predict good unlike LJ energy paramters for such mixtures. The reason 
is that all combining rules, except LB and M2, yield smaller values than the geometric mean. 
Such a lower unlike dispersion energy is supported by experimental noble gas mixture data 
[|2l |46l |571 158]|. However, for all other mixture types 2 to 5, the combining rules can predict 
better or worse than the simple LB combining rule. That all combining rules predict similar for 
mixture type 6 is due to the very similar values of both LJ parameters of the pure components. 
All in all, the simple LB combining rule has a RMSg of 5 % which is not significantly worse 
than those of the best rule (HMC) with 3.9 %, where a variation of ^ by 1 % influences the 
vapor pressure by about 3 %. 




(26) 
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Apart from noble gas mixtures, it can be concluded that none of the inverstigated combin- 
ing rules has a significant advantage over LB, which yields good bubble densities and vapor 
compositions, but vapor pressures that deviate from the experiment by an estimated 15 % on 
average. If highly accurate results for the vapor pressure are needed, an adjustment of the unlike 
LJ energy parameter to at least one experimental binary vapor pressure data point is necessary 
in most cases. 

5 Conclusion 

The dependence of vapor-liquid equilibrium properties on unlike LJ parameters was studied sys- 
tematically for the mixtures CO-1-C2H6 and N2-1-C3H6. This case study shows that the mixture 
bubble density is accurately obtained using the arithmetic mean of the like LJ size parameters 
as proposed by the Lorentz combining rule. The bubble density is insensitive to variations of 
the unlike LJ energy parameter. The vapor pressure is found to be dependent on both unlike 
LJ parameters. The same is found for vapor phase composition, but with a considerably lower 
sensitivity. Therefore, it can be recommended to use the Lorentz rule for the unlike LJ size 
parameter and to adjust the unlike LJ energy parameter to the vapor pressure. Recent investiga- 
tions on mixtures consisting of more complex molecules, such as ethanol Ii59l[60l , confirm that 
this procedure holds also in such cases. 

Eleven combining rules for unlike LJ parameters, some of which use also other properties 
than the pure component LJ parameters like ionization potentials or polarizabilities, were inves- 
tigated regarding their predictive power for vapor-liquid equilibria. They are compared to 44 
adjusted mixture models from previous work lfT3l[T4l[T5l . where the unlike LJ energy param- 
eters were optimized to yield experimental vapor pressure. It is found that the unlike LJ size 
parameters from the combining rules differ only very litte from the Lorentz rule. Many combin- 
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ing rules investigated here were fitted to noble gas mixture data and they indeed predict good 
unlike LJ energy paramters for such mixtures. However, none of the investigated combining 
rules has a significant advantage over LB for other mixture types. 
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Table 1 : Eleven selected combining rules together with the acronyms used here and correspond- 
ing references. In the last column, it is indicated whether the unlike LJ size parameter is given 
by the arithmetic mean, i.e. the Lorentz rule. 



Name 


Acronym 


Reference 


Lorentz 


Lorentz-Berthelot 


LB 


[22, 23] 


yes 


Kohler 


K 


126J 


yes 


Hudson-McCourbrey 


HMC 




yes 


Fender-Halsey 


FH 


m 


yes 


Hiza 


H 


[ISllllolEll 


no 


Sikora 


S 


m 


no 


Smith-Kong 


SK 


mm 


no 


Halgren 


HHG 


m 


no 


Waldman-Hagler 


WH 


Evil 


no 


MATAR-1 


Ml 


m 


no 


MATAR-2 


M2 


[48J 


no 
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Table 2: Binary mixture types (1-6) classified according to the molecular models of the two 
components. 





ICLJ 


2CLJQ 


2CLJD 


ICLJ 


1 


2 


3 


2CLJQ 




4 


5 


2CLJD 






6 
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Table 3: Unlike LJ energy parameters of the adjusted mixture models taken from lfT3l[T4l[T5l . 
grouped according to the mixture type as defined by Tabled 



Mixture 


Type 


Kb/kB]/K 


Mixture 


Type 


r adj / , 1 / 


Ne + Ar 


1 


53.752 


CO2 + C2H6 


4 


128.878 


Ar + Kr 


1 


136.280 


CO2 + CS2 


4 


170.086 


Ar +CH4 


1 


126.974 


CO2 + CI2 


4 


137.020 


Kr +Xe 


1 


190.225 


C2H4 +C2H6 


4 


106.470 


Kr +CH4 


1 


156.184 


C2H4 + C2H2 


4 


76.446 


Ne +N2 


2 


32.513 


C2H6 + C2H2 


4 


101.267 


Ne +O2 


2 


35.249 


C2H4 +CO2 


4 


95.579 


Ne +CO2 


2 


75.559 


C2F6 +CO2 


4 


105.045 


Ar +N2 


2 


64.065 


CO2 + Propylene 


4 


130.390 


Ar +O2 


2 


70.164 


C2H6 + Propylene 


4 


146.450 


Ar +CO2 


2 


124.610 


C2H4 + Propylene 


4 


109.654 


Kr +O2 


2 


82.030 


N2 + Propylene 


4 


69.491 


CH4 + N2 


2 


68.976 


C2H6 +C0 


5 


72.304 


CH4 + CO2 


2 


135.331 


CO2 + CO 


5 


75.719 


CH4 + C2H6 


2 


142.225 


N2 +C0 


5 


36.134 


CH4 + C2H4 


2 


109.268 


CS2 + R22 


5 


203.132 


CH4 + CO 


3 


74.256 


CO2 + R22 


5 


154.666 


N2 +O2 


4 


39.091 


CO2 + R23 


5 


127.914 


N2 +CO2 


4 


70.979 


R143a + R134a 


6 


168.985 


N2 +C2H6 


4 


67.344 


R125 +R143a 


6 


161.770 


N2 +C2H4 


4 


49.851 


R125 +R134a 


6 


168.663 


O2 +CO2 


4 


74.255 


R143a + R152a 


6 


177.304 
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Table 4: Deviations of the unlike LJ size parameters from the arithmetic mean, i.e. the Lorentz 
rule, for seven combining rules expressed by the root mean square RMS^ as defined by Equation 
(|25] ) for the different mixture types and for all 44 binary mixture models. 



Combining rule 


RMS^ / % 


Mixture type 


1 


2 


3 


4 


5 


6 


all 


S, Eq. (12) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


M2, Eq. dZU) 


0.3 


0.3 


0.2 


0.7 


0.1 


0.0 


0.5 


SK, Eq. dni) 


0.7 


0.6 


1.1 


0.5 


0.3 


0.0 


0.6 


HHG,Eq. ([H 


0.9 


1.0 


0.7 


1.7 


0.5 


0.0 


1.2 


H, Eq. d) 


1.0 


1.8 


0.3 


1.2 


0.5 


0.1 


1.3 


WH, Eq. dni) 


1.0 


1.2 


0.9 


2.0 


0.7 


0.0 


1.4 


Ml, Eq. (O 


1.0 


1.1 


0.8 


2.0 


0.6 


0.0 


1.4 
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Table 5: Deviations of the unlike LJ energy parameters from adjusted values for eleven combin- 
ing rules expressed by the root mean square RMSg as defined by Equation (|26l ) for the different 
mixture types and for all 44 binary mixture models. In the last column, it is indicated whether 
the unlike LJ size parameter is given by the arithmetic mean, i.e. the Lorentz rule. 



Combining rule 


RMS^ / % 


Lorentz 


Mixture type 


1 


2 


3 


4 


5 


6 


all 


HMC, Eq. © 


4.9 


4.9 


1.6 


3.1 


4.3 


1.4 


3.9 


yes 


M2, Eq. (I22l) 


3.2 


5.2 


4.0 


3.3 


5.3 


1.4 


4.1 


no 


SK, Eq. ([H) 


0.8 


6.1 


7.8 


2.5 


4.7 


1.3 


4.1 


no 


S, Eq. W 


0.8 


7.4 


8.1 


3.7 


4.9 


1.4 


4.9 


no 


LB, Eq. © 


6.8 


5.1 


0.3 


5.4 


3.9 


1.4 


5.0 


yes 


K, Eq. @ 


4.4 


11.0 


1.7 


5.0 


4.6 


1.5 


6.7 


yes 


Ml, Eq. (EQl) 


0.7 


7.2 


6.2 


8.2 


6.7 


1.4 


6.8 


no 


WH, Eq. ([H 


0.8 


7.4 


6.6 


8.0 


7.0 


1.4 


6.8 


no 


HHG, Eq. ([HI) 


2.9 


9.1 


11.5 


6.8 


8.8 


1.4 


7.3 


no 


H, Eq. ® 


2.4 


13.6 


2.8 


7.3 


4.8 


2.0 


8.4 


no 


FH, Eq. © 


1.6 


13.9 


20.5 


10.0 


13.3 


0.7 


11.0 


yes 
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List of Figures 



1 Vapor pressure of CO+C2H6 at 223 K and bubble point carbon monoxide mole 
fraction of 0.1980 mol/mol: deviations of the simulation results to the experi- 
mental vapor pressure of 5.614 MPa. The surface is a fit to 25 simulation results 
with varying r\ and ^. Solid slim line: r\ and ^ combinations yielding the exper- 
imental value. Solid bold lines: r\ and ^ combinations deviating by ±3 % from 

the experiment. Lines in the r|-^ plane: projections of the lines from the surface. [30 

2 Bubble density of CO-1-C2H6 at 223 K and bubble point carbon monoxide mole 
fraction of 0.1980 mol/mol: deviations simulation results to the experimental 
bubble density of 16.499 mol/1. The surface is a fit to 25 simulation results with 
varying T| and ^. Solid slim line: r\ and ^ combinations yielding the experimen- 
tal value. Solid bold lines: r\ and ^ combinations deviating by ±0.5 % from the 
experiment. Lines in the T]-^ plane: projections of the lines from the surface. . [31 

3 Dew point carbon monoxide mole fraction of CO-1-C2H6 at 223 K and bubble 
point carbon monoxide mole fraction of 0.1980 mol/mol: deviations of simu- 
lation results to the experimental dew point carbon monoxide composition of 
0.8065 mol/mol. The surface is a fit to 25 simulation results with varying r\ and 
^. Solid slim line: r\ and ^ combinations yielding the experimental value. Solid 
bold line: r\ and ^ combinations deviating by +4 % from the experiment. Lines 

in the T]-^ plane: projections of the lines from the surface |32 
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Superposition of the projections from Figures [I]to[3](CO+C2H6). The thin lines 
indicate the combinations of r\ and ^ where simulation and experiment coincide 
for the different properties. The thick lines represent the assumed uncertainties 
of dp = ±3 %, 5p = ±0.5 % and dy^.^^ = +4 %. The predictions of the com- 
binining rules are given by: O LB, ■ K, A HMC, V H, T S, □ SK, O WH, ♦ Ml, 
+ M2. The remaining two combining rules are not within the scale. Shaded: 
target area of r\ and ^ combinations yielding experimental values within their 
uncertainties 

Plot analogous to Figure |4] for N2+C3H6 at 290 K and bubble point nitrogen 
mole fraction of 0.1146 mol/mol. The thin lines indicate the combinations of 
r\ and ^ where simulation and experiment coincide for the different properties. 
The thick lines represent the assumed uncertainties of 5;?=±3 %, 5p = ±1 % 
and dy^2 = +4 %. The predictions of the combinining rules are given by: O LB, 
■ K, A HMC, □ SK, O WH, ♦ Ml, + M2. The remaining four combining rules 
are not within the scale. Shaded: target area of r\ and ^ combinations yielding 
experimental values within their uncertainties 

Predicted unlike LJ energy parameters from eleven combining rules com- 
pared to adjusted values for the mixture types ICLJ+ICLJ (1), 1CLJ+2CLJQ 
(2) and 1CLJ+2CLJD (3): O LB, ■ K, A HMC, ▲ FH, V H, ▼ S, □ SK, • 
HHG,OWH, ♦Ml, + M2 
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Figure 1 : Schnabel et al. 
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Figure 2: Schnabel et al. 
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Figure 3: Schnabel et al. 
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Figure 4: Schnabel et al. 
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Figure 5: Schnabel et al. 
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Figure 6: Schnabel et al. 
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